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Abstract 



We extract pole positions for the 7V(1535) and ./V(1650) resonances using two 
different models. The positions are determined from fits to different subsets 
of the existing irN — > irN, ttN —* nN and 7p — > np data and found to be 
1515(10)-i85(15)MeV and 1660(10)-i65(10)MeV, when the data is described 
in terms of two poles. Sensitivity to the choice of fitted data is explored. 
The corresponding irir and nn residues of these poles are also extracted. 
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I. INTRODUCTION 

Properties of the iV(1535) are difficult to extract from ttN — > nN and 'yN — > 71N due to 
the nearby rjN threshold As a result, a number of recent analyses have been based on 
data from ir~p — > r^n and 7^ — > ryp. These studies and coupled-channel analyses including 
pion production data, generally find values for the iV(1535) pole position, mass, width and 
photo-decay amplitudes which differ from those obtained from pion production data alone 
M - J/J. While these more recent studies suggest that some iV(1535) properties should 
be revised, the modification of any single quantity is complicated due to correlations. An 
example is the £ p parameter used by Mukhopadhyay and collaborators ||. This combination 
of the photo-decay amplitude (A 1 / 2 ), total (T T ) and r/N (T^) widths is relatively stable, even 
though values of A1/2 and IV vary by factors of two. Manley || has also noted that near- 
threshold ir~p — > rjn data provide little sensitivity to different parameter choices. 



In Ref. |T(| a two channel if-matrix model was presented for S-wave tiN and r/N scat- 



tering up to a center-of-mass energy of about 1700 MeV. There the main motivation was 
to extract the eta-nucleon scattering length (a) and effective range (r ) and to determine 



their uncertainties allowed by the existing ttN — > irN |]TTJ, nN — > r/N and 7p — > rjp 
H data. Below, this model will now be used to estimate the energies and residues of the 
S'-wave nucleon resonances iV(1535) and iV(1650) as complex poles of the T-matrix. Any 
problems with the iV(1535) may carry over to the nearby A r (1650) resonance, as the prop- 
erties of these two resonances are extracted from the same (Sn) itN partial wave and the 
same photoproduction multipole. 

In the model of Ref. [|Kj two poles corresponding to these resonances were included in 
the .fT-matrix, and their energies were tuned along with other parameters to give a fit to 
the data. However, in principle, the positions of these T- and i^-matrix poles can be quite 
different. Furthermore, it is the T-matrix poles that are of physical significance - hence 



their tabulation in Ref. ||12|| . A second reason for determining T-matrix pole positions is 
the greater variation of Breit-Wigner parameters within different parameterization schemes. 
For each pole, we have also extracted the corresponding residue. 

The present study differs from most of those carried out previously in that we have 
explored the effect of using different models and fitting different data sets. We have also 
considered, for nN elastic scattering, the effect of fitting the original experimental data 
rather than the amplitudes extracted from these data. In the next section, we compare the 
model used in Ref. |10| to that used in the VPI analyses. These two models have been 
utilized in our fits. In section III, we show our results and consider the factors responsible 
for differences in the extracted resonance parameters. 

II. FORMALISM 

The model of Green and Wycech is fully described in Ref. [pi]]. Here we repeat only the 
main elements, in order to facilitate comparisons with the VPI analyses. Both models are 
based on a 3-channel K-matrix formalism. In Ref. |I(J, a narrow energy range was chosen 
in order to justify the neglect of partial waves beyond / = 0. In the VPI fits, higher partial 
waves were included in fits which spanned a much wider energy range. However, these fits, 
while employing a multi-channel formalism, were not constrained by 77-production data. 



In the fits of Ref. [0], S'-wave scattering was considered in a system consisting of the 
two channels irN and r]N - here denoted simply by the indices 7r and 77. Then the .fT-matrix 
and the corresponding T-matrix , which are related by T = K + iKqT, can be written as 
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y JX^ JX m J y l-iq v A vv l-iq n A m J 

where q^^ are the center-of-mass momenta of the two mesons in the two channels tt, rj. The 
channel scattering lengths Aij are expressed in terms of the i^-matrix elements as 

A m = K m + iK^q^/ (1 - iq^K^). (2) 

As discussed in Ref. [[H]], these K-matrices are designed to account directly for several 
observed features of the experimental data such as the presence of two S"-wave ttN resonances 
and allow both to have a coupling to the two-pion channel. The latter channel is not treated 
explicitly, but introduced by reducing a three channel i^-matrix for ttN ,r]N and 7T7riV into 
the two channel form in Eq. [I]. The resultant .fT-matrices in this two channel model are then 
as follows: 



-fV^ — Eq _ E T El _ E "T 1 _ iqsK;j3 , ^TTT) ^TTT) T Eq _ E T t J^j^Jf^ 



K m - K m + ^ _ E + 1 1 _ ,^ K ^ , [6) 



wWp K - T3(0) , 73(1) K _ ^(0)73(0) _,_ y^r (1)73(1) 

wnere A 33 - £ + -Er±h , -"^3 



E -E 1 Ei-E ' ™ E Q -E 1 Ex-E 



K 3 - i2z!^ (4) 

and g 3 is a three-body 7r7riV phase space. In all there were 9 parameters in the i^-matrices 
and one parameter for normalizing the photoproduction data. 

In the second (VPI) approach, a Chew-Mandlestam i^-matrix has been used [TI|| to 
couple the elastic tcN channel to two inelastic channels, rjN and 7rA (in an 1=2 state). One 
starts with a 3x3 matrix: 



K = K TV K m . (5) 




Following the methods outlined in Ref. [IT], the T-matrix is written in the form 

T = p l ' 2 K{l-CK)- l p l '\ (6) 
and abbreviated as T = p 1//2 Tp 1//2 . In this notation, the elastic T-matrix is given by 

where 
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K = K + CvnK ^ + c ** k Za ^ 
1 — C^nK^ 1 — C^aKaa 



Ci being a dispersion integral |TT[ of phase space factors over the appropriate unitarity cut, 



and p = ImC. Inelastic channels are given by 



(■\ _i_ n t 

— V ^7r7T - 1 7T.. , 

T ™ = (1-C,tf«) (9) 

III. FITS TO DATA AND AMPLITUDES 
In Ref. [[10] the 10 parameters were determined by fitting the ^-production data of 
Refs. and the energy dependent Sll ttN — > nN amplitudes of [|TTJ over the center- 



of-mass energy range 1350 < E cm < 1700 MeV. However, a better approach is to fit the 
ttN — > ttN experimental data directly, thus avoiding the intermediate step of extracting par- 
tial wave amplitudes. Since the above K- matrix formalism is designed only for T 7r7r (S'll), 



the other partial waves are in the form advocated in [[TTJ. The procedure is, therefore, to 
first fit with this latter form all of the ttN — > ttN data over the full energy range (2.1 GeV) 
utilized by the VPI analyses. This fit is referred to as solution VPI. Data are then refitted, 
using the form of Ref. |10| , over the energy range 1350 < E cm < 1700 MeV, along with the 



rj— production data, with the non-Sll amplitudes kept fixed. In this case only the parameters 
of the above i^-matrix model are adjusted. These fits are referred to as solutions GW1X, 
where X/100 denotes the parameter combination (7373(0) related to the 7T7riV" branching for 
the iV(1535). The value of X/100 was varied from 0.00 to 0.04, thus generating solutions 
GW10 to GW14. The results of these fits are given in Table |. 

We also considered the effect of modifying the form used in fitting the rj photoproduction 
data. As a first step, an additional energy dependence was added. This amounted to 
replacing A(phot) in Ref. Q by A(phot) + B(phot)[E c , m , - 1485]/100. However, this had 
little overall effect with B(phot) being an order of magnitude smaller than A(phot). A 
second two-parameter form 

A oc a (1 + iT vN ) + — T vN , (10) 



analogous to that used in pion photoproduction fll3| , was also used. In the above, a and (3 
were taken simply as constants. This form was labeled GW2X, with X retaining its earlier 
meaning, and was used to generate the results presented in Table [D]. The actual values 
for the 9 parameters are given in Table |TTT| for GW11 and GW21 - the solutions with the 



smallest \ ■ m Table [TV| these parameters are converted into the more conventional form 



of Ref. [ 1 1 . Here it is seen that for GW21 some of these parameters are very different 



from their on-energy-shell counterparts, whereas those for GW11 are very similar to the on- 



energy-shell parameters in Table I of Ref. [|Hj] . The errors quoted in this table from Ref. [10 
will be used later, when error estimates on the pole positions and residues are made. 

In order to find the poles Ep — iTp/2 of the T- matrix in Eq. [L], the energy E appearing 
in Eqs. § and || and in the momenta qv, q v and ^3 was everywhere converted into E — iT/2. 
It is a built-in feature of the present i^-matrix formalism that the poles are at the same 
positions in all three matrix elements T V7T , T m , T^. This has been checked and found to be 
so within 10 keV. 
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The results are given in Table [V| and compared with the current values in Ref . [ 12[| . There 



it is seen that our results are consistent with previous values - especially those of Ref. |14 
As could be expected, our error bars are smaller due to the improvement in, and quantity 
of, the experimental data now being analysed. From Table [V| it is also seen that, although 
there is a dependence of the pole positions on the 2-pion branching, the differences - for 
the range of branchings considered - are essentially covered by the statistical errors on the 
positions. 

In Table |VI] a corresponding comparison has been made for the moduli and phases of 
the residues of the T n7T poles. This table also shows the moduli and phases for the two T m 
poles. Again as a consistency check we confirm that the residues at the T nr/ poles are simply 
the square root of the T nn and T m residues. 

In addition to the above poles there is the possibility of having poles on other Riemann 
sheets - far from the physical region - that can be probed by systematically reversing the 
signs of q-x and q v . These addition poles are quite symmetric - a point that can be understood 
in the limit where each Kij is a single pole A /7?77/ (E — E). In this case the T-matrix reduces 
to T oc [E - E - i7 ff g w - i^ n q n \~ x . 

IV. CONCLUSIONS 

In this paper we have extracted pole positions for the iV(i535) and iV(i650) resonances 
using two different models - the results being given in Table 0. It is seen that the iV(1535) 
pole positions predicted by these two models agree within about 15 MeV, whereas some 
of the predictions of the earlier models ||14|| , and |16| are considerably different. The 
iV(1650) pole values cannot be directly compared, as the most recent VPI fits have further 



poles and zeroes. However, if one compares to the 2 Sll resonance fit of Ref. |18[], agreement 
with our iV(i650) values is much improved. The reasons for differences can be manifold: a) 
the models used in the analysis are different, b) different subsets of partial wave amplitudes 
are fitted, c) data versus amplitudes are fitted, d) only certain data sets are fitted e.g. only 
tcN — » ttN or only nN — > tcN plus tcN — > rjN etc., e) the energy ranges over which the data 
is fitted can differ. We explicitly considered one such possibility in our analysis, by including 
either the Sll ttN partial-wave amplitude or the nN data. The N(1535) pole position was 
found to be quite sensitive to this choice, shifting about 50 MeV higher if the partial-wave 
amplitudes were fitted. This sensitivity was also seen in the associated residues. 

These various alternatives question the reliability of attempting to extrapolate into the 
complex energy plane the T-matrix from a model that only fits a limited selection of data 
over a limited range of energies on the real energy axis. In view of this, it would be desirable 
to have quantitative estimates of the errors expected on these pole positions. Unfortunately, 
for those fits involving directly all of the nN — > ttN, it is difficult to get a meaningful 



estimate of such errors. However, in the less ambitious approach of Ref. |L0|, where only 
the irN — > ttN Sll amplitudes - taking into account their error bars - were fitted, the 
Minuit minimization procedure gave error bars on the 10 parameters defining the model 
i.e. the Ei and ji in Eqs. ||,f|. Therefore, the errors 5Ep and 5Tp/2 on the pole positions 
Ep — iTp/2 could be obtained by repeating the calculation a large number of times for a 
random selection of the nine parameters defining the K- matrices of eq. ([l]) - as discussed in 
Ref. [|10j. This results in: 

5£ P (1535) « 10 MeV, <5r P /2(1535) « 10 MeV, <5£ P (1650) « 5 MeV, 5r P /2(1650) « 5 
MeV - values that were not very dependent on the actual pole positions. Such estimates 
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for SEp, 5Tp/2 are consistent with the spread of Ep — iTp/2 values from the various fits 
GW1X and GW2X. They are also very close to the correlated error estimates listed in Table 
V. Furthermore, it is seen that the position of the first pole, as given in the VPI and GW 
models, is consistent within these errors. Therefore, if we were to quote a single "best" 
number for the pole positions involving only two poles, Tables [VJ and [VJ] suggest: 
£p-TV/2(1535) = 1510(10) — «'85(15) and E P -TV/2(1650) = 1660(10) -i70(10) and the 
corresponding residues [\r\,9] for T n7T being [50(10), 0(10)] and [45(10), -40(10)]. However, 
the residues for T m depend strongly on the fit with the components for GW11 and GW21 
differing by about a factor of two and with the VPI estimate being somewhat closer to that 
ofGW21. 

In the above analysis the question of uniqueness arises. In the first model, the forms 
of Eqs. ^ and |4] are chosen with the physical idea in mind that there should be two basic 
resonances, which are compact in space ( as in a quark model) and so may be expected to 
be well represented by a pole in the K-matrix with a constant residue. Less compact objects 
would then need a form factor in place of the constant residue. This inclusion of explicit 
poles in the i^-matrix essentially guarantees poles in the T-matrix in the vicinity of those in 
the .ff -matrix. In the second of our models, poles in the i^-matrix can arise as a dynamical 
effect through coupling to high lying closed channels as in Eq. [5]. This alternative has also 
been discussed in Ref. fI7| , where the iV(1535) is treated as a KA bound state. This type 
of ambiguity has a long history and has been discussed in most detail for the interpretation 
of the A(1405) - see Refs. [fL9| . However, as emphasized in Ref. p0 |, the truth is probably 



somewhere in between the two above possibilities, with both mechanisms playing a role. 
This seems to be supported in Ref. |2~I||, where the authors conclude that the iV(1535) is not 
only generated by coupling to higher channels but " appears to require a genuine three-quark 
component" . In principle, with perfect data in all the relevant channels, the T-matrix should 
be highly constrained, so that only one prescription would succeed. However, in practice, 
the data have error bars and only cover a limited range, so that both approaches could give 
a fit to some of the data, but yield different pole positions. As a next step in resolving this 
uncertainty, all of the available data in ttN — > tcN, ttN — > r]N and rjN — > rjN (from final 
state interaction data in, for example, 7p — > rjp) should be treated simultaneously and not 
simply the selections used above. 

The authors thank Professors G. Ho filer for useful correspondence and also one of us 
(S.W.) wishes to acknowledge the hospitality of the Helsinki Institute of Physics. The 
authors were saddened to hear of the untimely death of Mijo Batinic, who has played an 
active role in the field of 77-nucleon physics. 
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TABLES 



TABLE I. Dependence of the fit, using the form of Ref. |T(]], with variations of the branching 
to ttttN. Notation is GW1X, where the parameter combination (see text) (/373(C)) takes on the 
value X/100. x 2 values are given for the fitted (2771) vr^p, (452) charge-exchange (CXS), (53) r\ 
photoproduction (Krusche), (11) -K~p — > rjn total cross section (Nefkens) data. The Sll column 
shows how well this fit to data reproduces the (60) VPI Sll single-energy points. The 2tt column 
shows the corresponding 2-pion branching ratio as a percentage. 



Soln. 


Total 


7T p 


CXS 


Krusche 


Nefkens 


Sll 2vr 


GW10 


7719 


6153 


1447 


66 


51 


105 


GW11 


7671 


6169 


1410 


62 


30 


98 2.6 


GW12 


7718 


6245 


1393 


57 


22 


101 5.2 


GW13 


7784 


6323 


1377 


67 


16 


100 7.8 


GW14 


7861 


6401 


1374 


73 


13 


104 10.0 


TABLE II. Notation 


as in Table 


|. Here the two-term form (see text) has been used to 


fit i] photoproduction dat 


a. The VPI solution is 


a fit to the elastic irN scattering database from 


threshold to 2.1 GeV (with only forward dispersion relation constraints), including 


; the ^-production 


data. 














Soln. 


Total 


7T p 


CXS 


Krusche 


Nefkens 


Sll 2vr 


GW20 


7688 


6144 


1442 


51 


50 


126 


GW21 


7600 


6102 


1420 


52 


25 


121 1.1 


GW22 


7627 


6147 


1410 


53 


17 


122 2.3 


GW23 


7690 


6227 


1398 


50 


15 


125 3.5 


GW24 


7774 


6309 


1395 


52 


18 


128 4.8 


VPI 


7540 


6040 


1397 


53 


49 


85 
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TABLE III. The optimised parameters defining the if-matrices for GW11 and GW21 in Tables 
| and 0. 





GW11 


GW21 




GW11 


GW21 


K 


0.1078 


-0.8336 


7.(0) 


0.0640 


0.1220 


K 


0.0157 


-0.1051 


7.(1) 


0.1071 


0.0913 


E (MeV) 


1538.5 


1582.5 


7r? 


0.2283 


0.6027 


Ei(MeV) 


1681.6 


1678.8 


73(0)^ (1535) (MeV) 


1.97 


1.97 








7 3 (l)<73(1650)(MeV) 


18.1 


16.5 



TABLE IV. 
as in Ref. p| 



The parameters in Table I II expressed in terms of widths and branching ratios 



GW11 



GW21 



T(Total)(MeV) 
r](br) 
ir{br) 
r(rotaZ,l)(MeV) 
n(br, 1) 



151.6 
0.576 
0.398 
150.4 
0.769 



354.4 
0.663 
0.326 
133.3 
0.758 



TABLE V. The Real and Imaginary parts of poles (Ep — iTp/2) in the complex energy plane 
compared with those quoted in [12]. The numbers in parentheses correspond to the 3rd Sll pole, 



for which some evidence was found in Ref. [15]. 



Reference £ P (1535)(MeV) r P (1535)/2(MeV) £ P (1650)(MeV) r P (1650)/2(MeV) 



Arndt 
Hoehler 
Cutkosky gj] 
This paper 
VPI 



vpko yj 

GW10 
GW11 
GW20 
GW21 



1501 
1487 
1510±50 

1510±3 

1499 
1510±8 
1514±9 
1502±3 
1509±3 



62 

130±40 

73±3 

55 
87±5 
90±6 
80±3 
82±4 



1673(1689) 
1670 
1640±20 

1666(1668) 
1657 
1662±3 
1658±4 
1667±2 
1663±2 



41(96) 

82 
75±15 

41(147) 
80 
70±5 
69±5 
60±4 
60±4 



S 



TABLE VI. The Moduli(|r|) and Phases(f9) of the residues of the two poles in both T n7r and 
T m compared wih those quoted in Ref. p2]] . Residues for the VPI 1650 MeV resonance are not 
included, as the VPI fit has an added pole in this region. 



T 


Reference 


\r\ (1535) (MeV) 


0(1535)(°) 


\r |(1650) (MeV) 


0(165O)( ' 


T 


Arndt || 


31 


-12 


22(72) 


29(-85) 




Hoehler @ 


— 


— 


39 


-37 




Cutkosky || 


120±40 


15±45 


60±10 


-77±25 




This paper 










T 

7T7T 


VPI 


40 
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GW10 


53±10 


-1±10 


54±5 


-43±5 




GW11 


57 


1 


54 


-48 




GW20 


43±5 


-10±5 


42±5 


-32±5 




GW21 


45 


-5 


42 


-37 


T 


VPI 


41 


-85 








GW10 


91±20 


-53±10 


8±5 


122±10 




GW11 


98 


-48 


11 


127 




GW20 


43±5 


-120±5 


6±10 


14±15 




GW21 


41 


-121 


8 


15 
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